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Discovery of Gamma-ray Pulsations from the Transitional Redback PSR 

J122T—4853 
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S. Fort 9,10 , F. Carnilo 11 , J. Deneva 12 , B. Bhattacharyya 5 , B. W. Stap 7 pers 5 , M. Kerr 13 

ABSTRACT 

The 1.69 ms spin period of PSR J1227—4853 was recently discovered in radio obser¬ 
vations of the low-mass X-ray binary XSS J12270—4859 following the announcement of 
a possible transition to a rotation-powered millisecond pulsar state, inferred from de¬ 
creases in optical, X-ray, and gamma-ray flux from the source. We report the detection 
of significant (5<r) gamma-ray pulsations after the transition, at the known spin period, 
using ~1 year of data from the Large Area Telescope on board the Fermi Gamma-ray 
Space Telescope. The gamma-ray light curve of PSR J1227—4853 can be fit by one 
broad peak, which occurs at nearly the same phase as the main peak in the 1.4 GHz 
radio profile. The partial alignment of light-curve peaks in different wavebands suggests 
that at least some of the radio emission may originate at high altitude in the pulsar 
magnetosphere, in extended regions co-located with the gamma-ray emission site. We 
folded the LAT data at the orbital period, both pre- and post-transition, but find no 
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evidence for significant modulation of the gamma-ray flux. Analysis of the gamma-ray 
flux over the mission suggests an approximate transition time of 2012 November 30. 
Continued study of the pulsed emission and monitoring of PSR J1227—4853, and other 
known redback systems, for subsequent flux changes will increase our knowledge of the 
pulsar emission mechanism and transitioning systems. 

Subject headings: pulsars: individual (J1227—4853)-binary-gamma rays: stars 


1. INTRODUCTION 


Millisecond pulsars (MSPs) are thought to be old neutron stars that have reach ed short ro¬ 
tati onal periods (< 10 ms) as the result of accretion from a binary companion (e.g., Alpar et a! 


19821 ). This hypothesis is supported b y the detection of X-ray millise cond pulsations from some 
low-mass X-ray binaries (LMXBs, e.g.. IWiinands fc van der Klislll99S). However, the ex istence of 
isolated MSPs such as the first MS P ever discovered, PS R B1937+21 ( Backer et al. 1982 1. requires 
a separate formation channel (e.g., Ilvanova et al . 20081 ) or some mechanism by which the binary 
system can either be disrupted o r the companion des troyed. With the discovery of the first “black 
widow” pulsar (PSR B1957+20, Fruchte r et al. 198a ) the theory that the binary companion could 
be totally or partially ablated by the energetic pulsar wind, once the MSP turned on as a radio 
pulsar, gained pop ulari ty. H owever, o ther black widows detected in the Galactic field before ~2008 
(e.g., J2051—0827, iSt a p pers et al. 1 996a Jbl) were not as energetic and their companions did not ap¬ 
pear to be losing mass at a significant rate, complicatin g the ablation hypothesis. Several theories 
have been put forth as to how the ablation proceed s (e.g.. lR,uderman et al.lll989l : lHarding fc Gaisser 


199fll : iLevinson fc Eichlerill99ll : iTakata et al.ll2012l ). but one probl em with this scen ario, in general, 
was the relative lack of black widows outside of globular clusters ([King et al.1 20013) • 


Black widow systems are characterized by relatively short orbital periods (few hours) and 
extremely low-mass (~ 0.02 M & ) companions while their cousins t he redbacks have more massive 
(~ 0.2 Mq) non-degenerate companions and show radio eclipses ([Roberts 20111 ). Since the launch 
of the Fermi Gamma-ray Space Telescope in 2008, radio astronomers have increased the number of 
known black widow and redback systems outside of globular clusters, and normal MSPs, v ia targeted 


obser vations of unassociated gamma-ray sources with pulsar-like characteristics (e.g., IRav et al 


2012 ). There are now at least 18 and 8 black widow and redback systems known, respectively, in 
the Galactic field compared to 3 and 1 before Fermi. These discoveries lend credence to the idea 
that isolated MSPs can be formed via ablation and, with a larger sample outside of the complicated 
environments of globular clusters, provide an opportunity to better understand this process. 


Recently, three redback systems have been observed to transition between rotation-powered 
MSP and accretion-powered LMXB states, implying that MSPs may undergo multiple transitions 
into and out of accreting states late in the recycling process and providing the opportunity to study 
MSP evolution in more detail. PSR J1824—24521, in the globular cluster M28, underwent a month- 
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long X-ray outburst, with accretion-powered X-ray pulsations, and was no longer d etectable as a 


rotat io n-powered radio pulsa r before dimming and again being seen as a radio MSP (jPapitto et al 


2013). Stappers et ah ( 20141 ) reported the disappearance of radio pulsations from PSR J1023+0038 
near the end o f June 2013, accompanied by a five-fold increase in gamma-ray flux above 0.1 GeV. 
Patruno et al. (2014J) reported X-ray and optical observations of this MSP and argued that there 
was evidence for the presence o f an accretio n disk. Accretion-powered X-ray pulsations at the 
spin period were later detected ( Archibald et al. 2014 ). The LMXB XSS J1227 0—4859 has been 


spatially associated w ith an unidentified gamma-ray source (2FGL J1 2 27.7— 4853, Ide Martino et al 


2010; Hill et ah 2011). By analogy with PSR J1023+0038, Hill et ah (2011) hypothesized that this 
might be a transitional system, but without concurrent spectroscopic obs ervat i ons it w as unc lear if 
accretion was still occurring. The situation was clarified in late 2012 when lBassa et ahl (|2013 .l2014l) 
reported a drastic decrease of the X-ray and optical flux of XSS J12270—4859. They suggested 
that the source may have transitioned to a rotation-powered pulsar state and noted a possible 
decrease in the gamma-ray flux of the associated point source, approximately coincident with 
the estimated transition time. Radio observations with the Giant Metrewave Radio Telescope 
discovered a pulsar with a 1.69 ms spin period, in a binary system with an orbital period of 6.9 
hours, consistent with the position of XSS J12270—4859, giving the sourc e the designatio n PSR 
J1227—4853 and firmly establishing that a state change had been observed ([Roy et ah 20151 . 120141 ). 
Once the spin and orbital periods were known, it was possible to detect X-ray pulsations in archival 
XMM- Newton observ ations before the state transition, but only during sub-luminous accretion disk 
states (JPapitto et ah 2015j). 


Using ~1 year of data from the Fe rmi Large Area T elescope (LAT. lAtwood et ahll2009h after 
the transition and the timing solution of lRov et ah (20151), we report the detection of pulsed gamma 
rays from PSR J1227—4853 with a significance of 5.0o\ We also characterize the spectral behavior 
of the source from the start of the mission until 2012 December and revisit searches f or variab ility at 
the orbital period, which is slightly different than the optical period reported bv lBassa et ahl (2014). 
Finally, we compare our results to previous studies of the gamma-ray emission from this system, 
discuss the energetics, and provide thoughts on emission models before and after the transition 
based on our observations. 


2. ANALYSIS AND RESULTS 

2.1. DATA SELECTION AND PREPARATION 

The LAT is a pair-production telescope sensitive to gamma rays with energies from 20 MeV 
to more than 300 GeV and observes the gamma-ray sky with 2.4 sr field of view. At 1 GeV, the 
LAT has a 68% containment radius of ~1° and a near on-axis effective area of 7000 cm 2 . For a 
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detailed description of the on-orbit performance of the LAT, see Ackermann et al. ( 2012b ' *1. 

We selected events from the P7REP LAT dats@ corresponding to the SOURCE class recorded 
between 2008 August 4 and 2014 December 1, with reconstructed directions within 15° of the 
position of PSR J1227—4853, energies from 0.1 to 100 GeV, and zenith angles < 100°. Good 
time intervals were then selected corresponding to when the instrument was in nominal science 
operations mode, the rocking angle of the spacecraft did not exceed 52° or the limb of the Earth 
did not infringe upon the region of interest, and the data were flagged as good. Analysis of LAT 
data was performed using the Fermi ScienceTool^f] v9r34p2. 


2.2. SPECTRAL ANALYSIS 


We first performed a binned maximum likelihood analysis on a 20° x 20° region centered on 
the source position using the P7REP_S0URCE_V15 instrument response functions (IR Fs) spanning 


t he e ntire data set. All sources from the third Fermi LAT source catalogs (3FGL, Acero et al 


2015) within 25° of PSR J1227—4853 were included in the model of the region and all spectral 


parameters of sources within 8° with average significances > 10cr over four years were left free. We 
also left free the normalization parameters of sources within 10° that were flagged as significantly 
variable, even if they did not pass the average significance cut. The spectral parameters of all 
other sources from 3FGL were kept fixed. We moved the position of the 3FGL sour c e asso ciated 
with PSR J1227—4853 (3FGL J 1227.9—4854) to the timing position from Rov et al. ( 2015 ). The 
Galactic diffuse emission was modeled using the glLiem-v05-revl .fits model while the isotropic 
diffuse emission and residual background of misclassified cosmic rays were jointly modeled using 
the isosource-v05.txt templat^E The normalization parameters of both diffuse components were 
left free. 


We modeled the spectrum of PSR J1227—4853 as both a single power law (Equation [lj) and 
an exponentially-cutoff power law (Equation [2]): 


dN 

~dE 



(1) 

(2) 


1 See also http://www.slac.Stanford.edu/exp/glast/groups/canda/lat_Performance.htm 

'http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/Pass7REP_usage.html 

3 The Fermi ScienceTools can be downloaded at http://fermi.gsfc.nasa.gov/ssc/data/analysis/software/ 

4 See http://fermi.gsfc.nasa.gov/ssc/data/access/lat/4yr_catalog/ 

5 The diffuse models are available for download at http://fermi.gsfc.nasa.gov/ssc/data/ 
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Nq is the normalization, Eq = 0.445 GeV is the pivot energy from 3FGL, T is the photon index, and 
E c is the cutoff energy. The best-fit spectral parameters are given in column 2 of Tabled) We detect 


a point source at the position of PSR J1227—4853 with a likelihood test statistic ([M attox et al. 
199 61) TS = 1967 (for one degree of freedom, the significance of a point source is Vts). 


Using the likelihood ratio test, a single power-law shape is ruled out, in favor of an exponentially- 
cutoff power l aw, w ith a confidence level of ~5.5<r, as indicated by the value of TS cu t (defined as 


m 


Abdo et ah 2013|). Tabled) also reports the photon (Tioo) and energy (Gioo) fluxes integrated 


from 0.1 to 100 GeV. 


We then analyzed ~4.25 years of data up to 2012 November 12 (start of the estimated transition 
time window) and ~2 years after 2012 December 21 (end of the estimated transition window) 
separately. The results from the pre- and post-transition periods are given in columns 3 and 4 of 
Tabled) respectively. The photon flux above 100 MeV dropped by factor of ~3 after the transition, 
the low-energy photon index hardened significantly, and the cutoff energy decreased. 


For both the pre- and post-transition data sets, we calculated the spectral points in Figured] 
by performing binned likelihood fits in the energy bands shown. For each energy band, we started 
from the corresponding best-fit model and only freed the normalizations of sources within 6° of 
PSR J1227—4853 and with TS > 100 when fitting the full energy range. We modeled the spectrum 
of PSR J1227—4853 as a power law (Equation |T]) with fixed photon index T = 2. For each energy 
band, we required that PSR J1227—4853 be detected with TS > 4 (~ 2cr) and with a predicted 
number of counts > 4, else a 95% confidence-level upper limit is plotted instead. For each data set, 
we found the highest-energy event within the 95% containment radius (as a function of event energy 
and angle with respect to the LAT boresight, 6, and on the location of conversion in the LAT) 
and only fit up to the energy bin containing this energy, resulting in the pre-transition spectral 
measurements extending to higher energies. 


Takata et al.l (|2014f ) have proposed a model to explain the gamma-ray flux increase in the 


transition of MSP J1023+0038 (jStappers et al. 12014) )) in which the rotation-powered pulsar mech¬ 
anism remains active (see Section 13.21 for further discussion) and the increase in gamma-ray flux 
is due to the presence of a new, dominant component. In order to investigate this possibility for 
PSR J1227—4853 during the LMXB phase, we performed an additional binned maximum likelihood 
analysis of the pre-transition data set in which we had two point sources at the position of PSR 
J1227—4853. The spectrum of the first source was set to the best-fit, exponentially-cutoff, power- 
law model from the post-transition data set with only the normalization parameter left free. The 
second source was fit well by a simple power law with a photon index of 2.45 ± 0.08; the addition 
of an exponential cutoff was not statistically justified. There is no significant preference for the 
two-source fit over one source with a curved spectrum. In the two-source fit, the normalization 
parameter of the first source did not change significantly relative to that found from fitting the 
post-transition data set. 
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Table 1. Spectral Fit Results 


Parameter 

Full Data Set 

Pre-transition 

Post-transition 

N 0 (KT 11 cm~ 2 s" 1 MeV" 1 ) 

3.19±0.12+°;J 8 

3.78±0.16loJo 

1.86±0.23±{j;J£ 

r 

2.13±0.051q[o6 

2.18±0.06^Q;gg 

1.66±0.22±0.06 

E c (GeV) 

6.8±1.6±J'? 

7.5±2.2lS;| 

2.8±1.1±0.2 

F] 00 (10 -8 cnr 2 s _1 ) 

6.41±0.34t°;*j 

7.89±0.42l°;g 

2.61±0.52+°;^ 

Gioo (1CT 11 erg cnr 2 s _1 ) 

3.50±0.12l°;^ 

4.16±0.15l°;| 

1.86±0.19±°;J° 

TS 

1967 

1745 

266 

TS cu t 

31 

19 

19 


Note. — For all parameters, the first uncertainties are statistical, 1 a, and the 
second reflect systematic uncertainties in the LAT effective area as detailed in the 
text. 
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Fig. 1.— Gamma-ray spectra of PSR J1227—4853 pre-transition (black lines and square points) 
and post-transition (blue lines and triangular points). The dashed lines show the lcr confidence 
region around the best-fit model for the pre- and post-transition data. 


2.2.1. SYSTEMATIC UNCERTAINTIES 


The systematic uncertainties in the LAT effective area (A e ff) have been estimated to be 10% 
for log 10 (£/l MeV) < 2, 5% for log 10 (£/l MeV) = 2.75, and 1 0% for log in (£?/ l MeV ) > 4 
with li near extrap o lation , in log space, between these energies (lAckermann et al.1 l2012bl ). Fol¬ 
lowing |Abdo et ah (2013), we estimated the effects of these uncertainties on our derived spectral 
parameters by generating bracketing IRFs using a modified A e g given by, 


A hrack (E,9) = A eS (E,9)(l + err(E)B(E)), (3) 

where err(E) represents the systematic uncertainties and we used B{E ) = ±1 as the bracketing 
function to estimate systematic uncertainties on IVo and B(E ) = ± tanh(log 10 (Fl/T , o)/0.13) for 
















































r and Eq, again with Eg = 0.445 GeV. For each fit with different bracketing IRFs, we generated 
new sourcemap^l for point sources with free parameters but used the sourcemaps generated with 
the P7REP_SOURCE_V15 IRFs for the fixed sources and the diffuse components. This last step is 
important as the diffuse components are tailored to the P7REP data and nominal IRFs, and the 
parameters from the fixed sources are from fits using the same IRFs. To estimate the systematic 
uncertainties on Fioo and Gioo, we recalculated these values from each bracketing IRF fit and took 
the maximum excursions from the nominal values as the systematic uncertainty. These estimated 
systematic uncertainties are given as the second uncertainties on spectral parameters in Table [B 


2.3. LONG TERM SPECTRAL BEHAVIOR 


Using a 15° radius selection we performed unbinned likelihood fits in 30 day time bins in order 
to characterize the gamma-ray flux change of PSR J1227—4853. For each time bin we started from 
the best-fit model of the entire data set and modeled the spectrum of PSR J1227—4853 as a simple 
power law (Equation [l|) with both Nq and T free. We kept the same point sources free as in the fit 
of the full data set, but only allowed the normalization parameters to vary. In our first attempt, 
the time bin spanning 55702.66 MJD to 55732.66 MJD was found to have a flux above 100 MeV 
twice as high as the value obtained from fitting the pre-transition data alone. No other 30-day time 
bins were found with similar flux values. To further investigate this apparent high flux point, we 
used LAT ScienceTool gttsmap to make a 3°x3° TS map, using unbinned likelihood, centered on 
the pulsar position in this 30-day time bin without the pulsar in the model. The pulsar position 
is well outside the 99% confidence-level contour of the TS map peak. Using the LAT ScienceTool 
gtfindsrc, we localized this emission to right ascension 12 h 26 m 9?6 and declination —49°35 , 24 ,/ 
(J2000) with a 95% confidence-level radius of rgs = 10 6 7 , approximately 42' from the position 
of PSR J1227—4853. Investigating this new source pos ition, we found the blazar candidate PMN 
J1225—4936 (CRATES J1225—4936, Healev et al. 2007 ) only 4 7 away, well within rg$. 


A likelihood fit of the data for this 30-day time bin with PSR J1227—4853 and PMN J1225—4936 
using power-law spectra in the model detected both sources significantly, with the flux split approx¬ 
imately evenly between the two, thus bringing the flux of PSR J1227—4853 in line with neighboring 
bins. If we include the candidate blazar in the model of the region and fit over the entire data 
set, we find the source with TS = 4, explaining why there is no corresponding source in the 3FGL 
catalog. A 30-day flux light curve analysis of the candidate blazar suggests that this is the only time 
bin in which the source is detected at > 3 a confidence level. The candidate blazar was not found to 
be significant near the estimated transition time window. We computed similar TS maps in time 
bins when the candidate blazar was not significantly detected, during the pre-transition period, and 
found that the peak of the TS map was always consistent with the position of PSR J1227—4853 


6 For an explanation of LAT sourcemaps generation, see the LAT data analysis tutorial at 

http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/binned_likelihood_tutorial.html#computeSourceMaps 
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within the 95% confidence-level contour. Therefore, we conclude that PMN J1225—4936 is likely a 
transient gamma-ray blazar with one flare, reaching a photon flux of (9.8 ± 4.1) X 10” 8 cm” 2 s” 1 
with T = 2.4 ± 0.2, during the reported time interval. 


CJ 

oo 

CD 
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15 — 

10 


V .}. ' j ^ 1 



Fig. 2.— The top panel shows Fioo in 30 day time bins. Points are only shown if PSR J1227—4853 
was found with TS > 12 (~ 2a for two degrees of freedom) and at least four predicted counts, 
else a 95% confidence-level upper limit is shown (red); error bars are lcr statistical only from the 
likelihood fits. The black dashed line shows the Fioo value from fitting the entire data set. The 
blue dashed line shows the Fioo value from fitting only the pre-transition data. The pink dashed 
line shows the Fioo value from fitting only the post-transition data. The bottom panel shows the 
best-fit r value for time bins where PSR J1227—4853 was significantly detected. In both panel s, 
the diagonal hatched region shows the estimated transition time window from Bassa et al. (2014). 
When maximizing the likelihood in each bin, a point source corresponding to the blazar candidate 
PMN J1225—4936 was included in the model of the region. 


We then recomputed the 30-day flux light curve of PSR J1227—4853 with PMN J1225—4936 
included in the model with a fixed power-law index of 2.4. Figure [2] shows Fioo versus time as well 
as the best-fit T for those bins in which we do not report a flux upper limit. We note that while the 
source was not found to be significantly variable during the first two years of the Fermi mission, it 
is flagged as variable in 3FGL; however, in Figure [2] the apparent variability is somewhat reduced 
compared to the result when the candidate blazar is not included in the model. In order to assess 
if the source should still be considered significantly variable once the candidate blazar is included 
in the model, we computed a spectral variability index for the pre-transition data similar to the 
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flux variability index used in the Fermi LAT catalogs ( Nolan et al.1 2 012 ). To do this, we repeated 
the spectral analysis in each 30-day time bin with the spectrum of PSR J1227—4853 fixed to the 
best-fit power-law model over the pre-transition data, recorded the negative log likelihood values 
(— In £ nu ii,i) and compared these to the values from the fits with the pulsar spectral parameters 
free (— In jC, t ) to construct the spectral variability index as, 


TS Var = — 2 (^[(—hiA) - (—ln£ null ,i 


(4) 


Assuming that the variability index is distributed as a x 2 with 50 degrees of freedom (52 months 
minus 2 since we let both the index and normalization be free in the fits) we need TS Var > 76.15 
to say the source is variable to better than 99% confidence. Using Equation |4] we find TS Var = 
95, compared to TS Var = 125 without the candidate blazar in the model, suggesting that the 
spectrum is still variable during the pre-transition period but the amount of variability is reduced. 
It is important to note that Fermi LAT catalogs only allow the source normalization to be free 
when computing the variability index, whereas we let both the index and normalization be free, so 
while the catalog analysis addresses whether the flux is variable our analysis addresses whether the 
entire spectrum is variable. 


The LAT upper lim it approximately halfway through the transition time window indicated by 
the X-ray observations (IBas sa et al.1120141 ) could provide a more precise estimate of the transition 
epoch. Examination of a flux light curve in two-day bins (not shown) bracketing the transition time 
window indicates a significant drop in TS near 2012 November 30 with predominantly TS > 4 (as 
high as 26) before and TS < 4 (with one bin having TS = 5) after. We have verified that the drop 
in TS after this date is not due to decreased exposure. However, the flux appears to settle gradually 
in the post-transition time period with the emission starting softer than the pre-transition emission 
and gradually hardening. This may suggest that the transition from LMXB to rotation-powered 
pulsar state progresses over a longer time period than when transitioning in t he other direct ion, as 
implied by the abrupt increase in gamma-ray flux seen for PSR J1023+0038 (jStap pers et ah 20141). 


2.4. PULSED GAMMA RAYS 


To test for ganmra-ray pulsations from PSR J1227—4853 we kept only events within a 3° 
radius of the radio position af ter the tr ansi tion and assigned a r otational p hase to each LAT event 
using the timing solution of Rov et al. ( 2015 1 the fermi plugin ( Rav et al. 2011 1 to the Tempo2@ 


software (Hobbs et al.l 200a). We tested for pulsations in the LAT data from the start of radio 


monitoring (MJD 56 7 07.98) to the end of the da ta set described in Section 12.11 Using the H test 
(|de Jager et al.1 Il989l : Ide Jager Sz Biischind |201Cj) together with weights for each photon derived 


'http://sourceforge,net/projects/tempo2/ 
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from the spectral analysis (as formulated by lKerdl201ll h we obtain a test statistic value H = 35.3, 
implying a chance probability of 7.3xl0~ 7 , or a 5 .Oct significance, indicating a significant detection 
of gamma-ray pulsations. The H test is a powerful test for rejecting the null hypothesis of a uniform 
distribution of phases, and including the photon weights, which represents the probability a photon 
originated from the pulsar instead of a different source, further increases its sensitivit y. For in¬ 


depth descriptions of the test and i ts applicabi li ty to gamma-ray p ulsar data see Ide .lager et al 


(|1989|) ; Ide .Tager fc Biischind (120101 ) ; iKerrl (1201 ll ) ; lAbdo et al.l (120131 ) . As a check on our detection, 
we assigned random phases, from a uniform distribution between 0 and 1, to our calculated spectral 
weights and recalculated the H statistic. In 1x10 s trials, the largest H-test value we found was 
28.5, indicating that the probability of obtaining an H-test value > 28.5 by chance was < 1 x 10~ 5 , 
supporting the significance of our detection. 

We investigated how far backward in time, toward the end of the estimated transition time 
window, we could extend our data set for the gamma-ray light curve. The pulsed significance 
continued to increase when adding data back to approximately MJD 56650 (^2 months before 
the start of the ephemeris validity range). Adding events before this date made the detection less 
significant, indicating that the timing model is not sufficient to extrapolate backward beyond this 
point. This is consistent with the fact that the radio model determined from the 2014 February 
throug h December t iming data does not extrapolate backward to the 2013 November data, as found 
by Rov et al. ( 2015 b 


We also explored the possibility of extending the pulsar ephemeris back in time before the radio 
observations using the extra year of LAT gamma-ray data after the transition. An extended timing 
solution could lead to better constraints on the system, especially the orbital parameters, owing 
to the longer baseline and would address the question of whether the system has been an active 
gamma-ray pulsar since the state transition. However, the relatively low signal-to-noise ratio of the 
gamma-ray pulsations from PSR J1227—4853 precludes derivation of a reliable timing solution with 
the LAT. In particular, exploring the parameter space in the neighborhood of the radio ephemeris 
resulted in many local maxima, but each with large parameter uncertainties and too low statistical 
significance to be confidently trusted. Given the potential insights to be gained from being able to 
extrapolate the timing solution back further, this is cert ainly worth revisiting with more data and 
the forthcoming, more sensitive Pass 8 LAT data (lAtwood et al. 2013). 


The > 100 MeV light curve of PSR J1227—4853 using all events within a 3° radius of the radio 
position and detected between 56 650 M JD and the end of our data set is shown in Figure [3j The 
Parkes 1.4 GHz radio profile from I Roy et al. (2015j) is over-plotted. The absolute phase alignment 
was done by using Tempo2 to define pulse phase 0 using the 1.4 GHz radio TOAs, and plotting 
the Parkes radio template with the LAT light curve aligned to that fiducial phase using the fermi 
plugin for Tempo2. 
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Pulse Phase 


Fig. 3.— Phase-aligned gamma-ray (histogram, 20 bins per rotation) and 1.4 GHz Parkes (red line) 
light curves of PSR J1227—4853, with two rotations shown for clarity; the radio profile has been 
renormalized and the zero level adjusted for ease of viewing and has arbitrary units. The statistical 
uncertainties for the LAT light curve bins and the gamma- ray backgro u nd le vel, horizontal dotted 
line, are derived from the spectral weights estimated as in lAbdo et alj ({2013) ■ The low-level peak 
at phase ~0.5 in the Parkes light curve is a real interpulse, which becomes dominant at lower 
frequencies. 


The gamma-ray light curve of PSR J1227—4853 can be satisfactorily fit with one broad Gaus¬ 
sian, FWHM = 0.38 ± 0.09, at 0.09 ± 0.03 in phase (assuming the estimated background level 
from the weights to reduce the number of fit parameters). A fit with two closely-spaced Gaussian 
peaks is also acceptable, but the improvement over the single-peak fit is not significant with the 
current statistics. Figure H] shows the energy evolution of the gamma-ray light curve. When only 
considering events with energies from 0.3 to 1.0 GeV, the light curve peak is almost at phase 0 
while when looking at the light curve for events with energies > 3 GeV the peak is sharper and 
near 0.1 in phase. This behavi or is si milar to w hat is seen for other pulsars with two closely-spaced 
peaks (e.g., PSR J0007+7303, Abdo et al. 2013 1. 
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Pulse Phase 

Fig. 4.— Gamma-ray light curves of PSR J1227—4853 (15 bins per rotation) in 4 energy bands as 
labeled, with two rotations shown for clarity. The horizontal dotted lines in each panel show the 
estimated background levels. Uncertainties and background levels are estimated as for Figure [3j 















































































































































































































































-14- 


T he gamma -ray peak is nearly aligned with the main peak in the 1.4 GHz Parkes radio profile. 


Rov et al.1 (2015)) also presented the pulse profile of PSR J1227—4853 at 607 and 322 MHz. The 
lower-level peak in the 1.4 GHz profile, near ^0.5 in phase, is a real feature and becomes dominant 
at lower frequencies. In fact, at 322 MHz the main 1.4 GHz peak disappears entirely. 


2.5. TEST FOR ORBITAL MODULATION 


In order to test for modulation at the orbital period, we first made a counts light curve in 30 
second time bins from the post-transition data set and then calculated the orbital phase and LAT 
exposure for each time bin. This allows us to correct for potential exposure variations across the 
orbital period, which has been show n to be a necessary step to avoid detection of false modulation 


(Kerr 


2010 


Ackermann et al.1 2012a ). 


We created a new null distribution (‘hnuii) for the variability tests by binning the exposure in 
1000 bins of orbital phase and normalizing. This 4> nu n reflects the variation with orbital phase we 
would expect for a non-varying source due to differences in exposure. Following Kerri (|2010l ). we 
used our new 4>nuii and the spectral weights calculated in Section 12.41 in a weighted H test and Z 2 m 
test with two harmonics to test for modulation at the orbital period. The significance values from 
the weighted H and tests were l.lcr and 1.4<r, respectively, indicating no strong evidence for 
modulation at the orbital period. 


We split the post-transition data into 10 orbital phase bins, using the binary period from the 
radio timing solution. The good time intervals for each bin were corrected to properly account for 
the orbital phase selections. Starting from the best-fit model with the spectrum of PSR J1227—4853 
modeled as an exponentially cutoff power law (Equation [2]), we kept only the normalization pa¬ 
rameters of sources free and performed binned likelihood fits in each bin. The resulting orbital flux 
light curve is shown in Figure [5j 

Using the phase-averaged flux in a y 2 analysis (neglecting the uncertainty on the phase- 
averaged value) gives a reduced y 2 value of 1.47 (with 9 degrees of freedom), a larger y 2 value 
would be found 15% of the time by chance, suggesting this is an acceptable fit. This reduced y 2 
value suggests that any variations in Figure[5]are consistent with statistical fluctuations, supporting 
the conclusion that the gamma-ray emission from PSR J1227—4853 is not modulated at the orbital 
period. 


3. DISCUSSION 


Hill et al. (201 lj; used ~2 years of LAT Pass 6 data to analyze the region around PSR 
J1227—4853 while it was still, we now recognize, in the LMXB phase. Our best-fit T and de¬ 
rived photon flux for the pre-transition data set agree with their results, but our best-fit cutoff 
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Fig. 5.— Flux light curve versus orbital phase for PSR J1227—4853 for the full post-transition 
data set; uncertainties are la statistical only from the likelihood fits. Each point was found with a 
significance > 2cr. The dashed black line is the phase-averaged flux given in Tableland the dashed 
red lines show ±1 a (statistical only). 


energy is nearly twice their value of Eq = 4.1 ± 1.3 GeV, though consistent within the joint 
uncertainties. 


Xing fc Wand (2014) present a similar LAT study of PSR J1227—4853 using ~6 years of 
data. Their spectral analysis was restricted to events with energies > 0.2 GeV. Using our best 


fit, we derive a flux from 0.2 to 100 GeV of (2.79 ± 0.96) x 10” 8 cm -2 s" 1 over the entire data 
set, statistical uncertainty only; this value and our best-fit T and Eq are in good agreement with 
their results. Additionally, our best-fit, pre-transition model gives a flux from 0.2 to 100 GeV of 
(3.34 ± 0.12) x 10~ 8 cm" 2 s -1 , and our best post-fit transition model gives (1.44 ± 0.17) x 10~ 8 
cm -2 s _1 , statistical uncertainties only. Both of these values and the other spectral parameters are 
in good agreement with their results. 


A direct comparison of our flux light curve to that of Xing &; Wang (201411 is complicated by 


their use of 0.2 GeV as a minimum energy and by the fact that they do not impose a cut on TS 
when deciding whether to plot points or upper limits in the flux light curve. The spectrum is still 
variable during the pre-transition data but the significance of the variability is reduced in our result 
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due to the i nclusion of t he s ource we associate with PMN J1225—4936. The > 0.2 GeV flux light 
curve of Xing Jz Wang (20141) does show a similar increase around the time of our claimed flare 
from this candidate blazar. 


Xing fc Wand (120141 ) have claimed marginal ev idence for m o dulat ion of the gamma-ray emis 


sion from PSR J122T—4853 at the optical period of lBassa et al. (2.Q14)) after the transition, which 
is not seen in our analysis. Using the optical period and 10 orbital phase bins, we found that the 
exposure varies by < 2% from the average in any given bin, indicating that their claimed orbital 
modulation is not likely due to varying exposure alone. We note, however, that their stated H 
test value of ~10 corresponds to a 2.4<r detection, not 3cr as claimed, and trials factors must be 
accounted for, reducing the significance even further. In attempting to verify this claim we folded 
the same time interval of LAT data used in their paper at the optical period and used the same 
energy and radius cuts, but could only recover an unweighted H-test value of 8.1 (2.1<r). Our orbital 
flux light curve (Figure [5]) does have two bins that are 1.5~2 cj below the phase-averaged flux value, 
but both are well away from superior conjunction (<^ or b = 0.25) and thus are not easily explained 
as eclipses of the gamma-ray em ission site by t he co mpanion (a proposed explanation for the dip 
seen in the orbital light curve oflXing Sz Wand 2014). 


3.1. ENERGETICS 


An important quantity for understanding the pulsar emission mechanism is the gamma-ray 
luminosity, defined as L 7 = AnfnGiood 2 , where fn is a beaming co rrection factor accounting 
for the non-isotropic nature of the emission (e.g., Watters et al.1 20091 ) and d is the distance to 
the pulsar. T he r adio timing solution gives a dispersion measure (DM) of 43.4 pc cm~ 3 for_PSR 
J122 7—4853 ( Rov et al. 2015 ).^Using the NE2001 electron-density distribution ( Cordes fc Lazio 
20021 ) yields d = 1.4 ± 0.2 


Johnson et al.l (|2014l ) fit the g amma-ray and ra dio light curves of all MSPs in the second 


LAT catalog of gamma-ray pulsars (|Abdo et ah 20131). using models that assumed the bulk of the 
gamma-ray emission comes from regions high in the pulsar magnetosphere. They found that, for 
configurations capable of producing light curves that match the observations, fn values are typically 
close to but slightly less than 1. For the remainder of the discussion we will assume fn = 1. 


Using the post-transitio n spectral results and the DM distance, we find L 7 = (4.4 ± 0.5 ± 
1.3) x 10 33 erg s _1 (following lAbdo et al. 20131. we separate the uncertainty due to the gamma-ray 


spectral fit, first value, and that due to the distance, second value). We define the gamma-ray 


s The distance uncertainty is calculated by assuming a ±20% uncertainty on the DM, which is much larger 
than the formal DM uncertainty from the timing solution and is meant to reflect uncertainties in the electron- 
density distribution model, and calculating the resulting change in d using the online NE2001 implementation 
http://www.nrl.navy.mil/rsd/R0RF/ne2001/ 
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efficiency as ?? 7 = (L^/E) x 10 0%, where E = 47r 2 TP/P 3 = 0.9 x 10 3j erg s _1 is the spin-down 
luminosity of PSR J1227—4853 (Roy et aL 2015), / is the moment of inertia and is taken to be 10 45 
g cm 2 , P is the spin period of the pulsar in seconds, and P is the time derivative of the pulsar spin 
period. Using this definition and the values for PSR J1227—4853 yields r/ 7 = (4.9 ± 0.6 ± 1.4)%. 
Comparing to other gamma-ray pulsars with similar E , an effi ciency of ^5% is on the low end but 
consistent with the L 7 vs. E trend seen for LAT pulsars (lAbdo et al. 2013). 


It is also of interest to examine the energetics of the pre-transition emission. Similar consider¬ 
ations give L 7 = (9.8 ± 0.4 ± 2.8) x 10 33 erg s” 1 and r/ 7 = (10.9 ± 0.4 ± 3.1)%, indicating that 
there is sufficient rotational energy to power that emission as well. 


3.2. EMISSION MODELS 


3.2.1. PULSED LIGHT CURVES 


The prevailing mo dels for high-energy pulsar gamma-ray emission , such as the outer gap (OG 


e.g.. iCheng et al.l Il986l) and slot gap (SG. iMuslimov Sz Hardina 120031 120041 ) models, assume that 
particles are accelerated in vacuum gaps along curved magnetic field lines and emit photons in the 
outer magnetosphere, far from the stellar surface. In the OG model, particles are assumed to be 
accelerated and emit gamma rays via curvature radiation in a vacuum gap above the null-charge 
surface (the geometric surface defined by the condition Q ■ B = 0), bordered by the last-closed 
field lines (those which close at the light cylinder and are able to rotate with the pulsar). The 
SG model, on the other hand, assumes particles are accelerated in gaps bordering the last-closed 
field lines, emitting gamma rays along the way, from the stellar surface to near the light cylinder. 
In both models, the bright, and typically sharp, peaks observed in gamma-ray pulsar light curves 
are caustics where time-of-flight delays and relativistic aberration result in photons emi tted at 
different alti tu des on trailing field lines arriving at an observer coincident in phase (e.g., Mprini 


1983 


Venter et al. 20091 ). For both the SG and OG models, the vacuum gaps are assumed to be 


relatively narrow, a few percent of the polar cap size, as the rest of the magnetosphere is thought 
to be filled with ch arges th at shor t out the accelerating electric field. The pair-starved polar cap 


model (PSPC, e.g., Hardin g et al. 2005 ). however, was developed for MSPs and old non-recycled 


pulsars where it was thought that the weaker magnetic fields would not be capable of producing 
enough charges, through pair cascades, to screen the accelerating field over the entire polar cap. 
This means that the entire open field line region above the polar cap is available to accelerate 
particles and produce gamma rays. 

With standard, low-altitude models of pulsar radio emission, the OG, SG, and PSPC models 
predict significant phase lags between features in the radio and gamma-ray light curves, observed for 


all n on-recycled pulsars (except the Crab) and most MSPs known to emit gamma rays (lAbdo et al 
2013 ). However, a subclass of gamma-ray MSPs has been established in which the gamma-ray 


and radio peaks occur at (nearly) the same rotational phase. One explanation for this subclass 
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has been that the radio emission is also produced far from the magnetosphere, co-located with the 
gam ma-ray emission region, described by al ti tude-limited version s of the the OG and SG models 
(e.g., IVenter et al.ll2012l : Ijohnson et al.ll20141 h IVenter et al.l ( 2012 ) also introduced the low-altitude 
slot gap (laSG) model for this subclass of gamma-ray MSPs. In the laSG model, the gamma rays 
and radio are both produced within a few stellar radii of the neutron star surface in a narrow 
gap near the l ast-cl osed field lines. Similar to PSR J1824—2452A (in the globular cluster M28, 
Johnson et ah 2013|) the radio and gamma-ray light curves of PSR J1227—4853 appear to have 
some aligned and some non-aligned features, so PSR J1227—4853 does not cleanly fit into this 
gamma-ray MSP subclass, though the short spin period and large E would make it, a priori, a 
prime candidate. 


Johnson et al.l (J2014I) jointly fit the radio and gamma-ray light curves of MSPs in the second 
LAT cat a log of gamma-ray pulsars using geometric realizations of the OG, two-pole caustic (TPC, 
Dvks fe Rudakll2003l . taken to be a geometric realizati on of the SG), PS PC, and laSG gamma-ray 
emission models and hollow-cone plus core beam (e.g., Story et ah 2007)), altitude-limited versions 
of the OG and TPC models, or laSG radio emission models. With the current statistics for PSR 
J1227—4853, a constraining fit is not possible, but qualitative comments can be made concerning 
the probable system geometry and viable emission models. The emission geometry of a pulsar is 
typically defined by two angles, the angle between the spin and magnetic axes (a) and the angle 
between the magnetic axis and the observer’s line of sight (/?). Often, it is instead desirable to work 
in terms of a and the angle between the spin axis and the observer’s line of sight (£ = ft — a). 
For pulsars in binary systems, another important angle is that between the orbital axis and the 
observer’s line of sight (*). 


de Martino et al. 1120141 1 used optical observations of XSS J12270—4859 to constrain the orbital 
inclination to be 45° < i < 65°. The recycling process is thought to align the orbit and spin axes, 
such that we can assume i ~ £, translating the previous constraints into limits on £. Similarly, 
t he o bserva tion of radio eclipses for significant portions of the orbit near superior conjunction 
(iRov et al. 20151 ) suggest i > 60°. The radio profile at 1.4 GHz (Figure [3]) consists of two peaks 
separated by approximately 0.5 in phase. Assuming a hollow-cone beam and £ consistent with the 
optical constraints on i, the radio profile suggests a near 90°, but not exactly 90° as the two peaks 
are not symmetric. However, standard outer-magnetospheric gamma-ray emission models cannot 
produce the observed peak separation and phase lag with respect to the radio for those geometries. 
Acceptable simulated gamma-ray light curves can be found using the OG, TPC, or PSPC models 
and a hollow-cone radio beam with reasonable values of £, but all require low values of a (< 20°) 
and predict only one radio peak. If we fit only the gamma-ray light curves, the TPC model finds a 
best-fit geometry with £ near 65°, agreeing well with constraints on i, and a near 30°, while the OG 
model prefers large a values near 80° and small £ near 15°, not in agreement with the constraints 
on i. We remind the reader, however, that with the current gamma-ray statistics it is difficult to 
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make definitive statements regarding emission geometries that can reproduce the observed light 
curves of PSR J1227—4853. 


Another possibility is that the 1.4 GHz radio peak that is nearly aligned with the gamma-ray 
peak originates in an extended region in the outer magnetosphere co-located with the gamma- 
ray emission region, while the other radio peak arises from lower altitude nearer the polar cap. 
This model would eliminate the need for a to be near 90° in order for two conventional radio 
peaks to be visible. The gamma-ray only TPC geometry with a near 30° and £ near 65° would 
then be consistent with the orbital inclination constraint and with single gamma-ray and radio 
caustic peaks. Different locations for the two radio peaks could also account for their very different 
spectral evolution. Radio polarization measurements could verify this model, which would predict 
depolarization and rapid position angle swings in the caustic peaks but higher levels of polarization 
in the other peak. 


Du et al. (120151 ) reproduced the 1.4 GHz radio, X-ray, and gamma-ray light curves of PSR 


J1824—2452A using a single-pole annular gap model. The emission regions for all three wavebands 
were modeled as extended in the pulsar magnetosphere with some overlap, explaining the fact that 
some, but not all, of the peaks are aligned in phase. A similar model may be able to explain 
the radio and gamma-ray light curves of PSR J1227—4853, though an independent constraint on 
a would be helpful for such an endeavor as those authors have not yet developed a statistical 
technique to choose the best-fit geometry from profile fits. 


3.2.2. PRE- TRANSITION EMISSION 


Takata et ah) (2014) analyzed multi-wavelength data for PSR J1023+0038, characterizing the 
transition of the system from rotation to accretion-powered state that occurred near the end of 2013 
June. They proposed an emission model in which an accretion disk forms but does not penetrate 
inside the light cylinder (defined by the cylindrical radius at which rotation with the star requires 
moving at the speed of light) suggesting that the rotation-powered pulsar emission mechanisms are 
still active. In their model, the observed gamma-ray flux increase represents the emergence of an 
additional component, not modulated at the orbital period, due to inverse Compton scattering of 
the pulsar wind off of ultraviolet photons from the accretion disk. From the two-source fit of the 
PSR J1227—4853 pre-transition data detailed in Section \2 .21 we cannot rule out that the rotation- 
powered pulsar mechanism was active while the source was in the LMXB state. However, without 
a detection of gamma-ray pulsations during this time period we cannot confidently say that it was 
active either. 


Papitto et al. (120151 ) argued that the X-ray pulsations at the spin period during the LMXB 


phase from PSR J1227—4853 are most plausibly explaine d as a ccretion rather than rotation- 
powered, which argues against the model of T akata et al, ((2014) to explain the pre-transition 
gamma-ray emission. Recent observations of PSR J1023+0038 during the LMXB state also ar- 
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gue against this model. Archibald et ah (120141 ) have detected coherent X-ray pulsations, most 
likely accretion powered, from PSR J1023+0038 after the recent state transition. Radio observa¬ 
tions, at frequencies up to ~5 GHz, h ave failed t o detect pulsations at the rotational period while 
the source i s in th e LMXB state (Bogdanov et alJ 12014 ). Finally, radio imaging observations by 
Deller et al. (2014), during the LMXB state, found a variable, flat-spectrum radio source that is 
incompatible with the previously very-steep spectrum of PSR J1023+0038. 


Papitto et ah (2014]) proposed a propellor model to explain the pre-transition (LMXB phase) 


X- and gamma-ray emission from PSR J1227—4853. In their model, the accretion disk penetrates 
inside the light cylinder (with an inner disk radius of ~40 km for PSR J1227—4853) quenching 
the rotation-powered pulsar emission mechanism. The observed X- and gamma-ray emission is 
the result of synchrotron and synchrotron self-Compton scattering from electrons accelerated at 
the interface between the disk and the puls ar magnetosphere . For reasonable parameters, their 
model was able to reproduce the spectrum of Hill et all (120111). Both the spin period and inferred 
dipole mag netic field stre ngth reported bv IRov et all (2015) are in good agreement with the values 
assumed bv lPapitto et al.1 (12014 ), su ggesting that their model is a viable explanation of the emission 
during the LMXB phase. [Papitto et al . (2015! ) inferred an accretion disk truncation radius for XSS 
J12270—4859, from the luminosity of the X-ray pulsations, that was outside of the corotation radius, 
which would prevent accretion. To reconcile this with the apparent accretion-powered nature of 
their observed X-ray pulsations, they theorized that the disk mass accretion rate was actually ~30 
times larger than their inferred value, but that only a small portion was channeled to the neutron 
star surface with the rest being eje cted i n an outflow, a scenario that is in agreement with the radio 
spectrum observed by iHill et al.1 ( 201 ll ). Such an outflow would argue in favor of the propellor 
model. 


4. CONCLUSIONS 


It is now clear that PSR J1227—4853 is one of three redback MSPs, to date, observed to have 
transitioned between an LMXB state and a rotation-powered pulsar. We have detected significant 
gamma-ray pulsations from PSR J1227—4853 after the transition in late 2012 and characterized the 
spectral differences pre- and post-transition. Analysis of the LAT data during the transition time 
window in 2-day time bins suggests that the transition occurred near 2012 November 30. After the 
transition the photon (energy) flux dropped by a factor of ^3 (~2), the photon index hardened, 
and the cutoff energy decreased. We find no evidence in the LAT data, pre- and post-transition, 
for significant modulation at the orbital period. PSR J1023+0038 has been observed to transition 
in the other direction, from rotation-powered pulsar to LMXB, displ aying a larger gamma -ray flux 
i ncrease, with estim ates of the change varying from a factor of ~5 (IStappers et al.l 120141 ) to ~10 
(ITakata et al.1 1 20141 ). PSR J1023+0038 also displays a lower cutoff energy and harder photon index 
during the rotation-powered state compared to the LMXB state ( Takata et al.1 120141 1. During the 
rotation-powered pulsar state, the gamma-ray flux from PSR J1227—4853 is times greater than 
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that from PSR < 11023+ 0038, which is at a similar distance (1.4 kpc) but has a lower E (4.4x 10 34 erg 


s 4 ) ( Archibald et ah 2013). While the viewing geometry is not precisely known for either system, 


we can speculate that the difference in pulsed fl ux is likely due to beaming and a less favorable 
geometry for PSR J1023+0038 (e.g., Fig. 20 of Venter et al. 2009 1. During the LMXB state the 
observed gamma-ray fluxes from PSR J1227—4853 and PSR J1023+0038 are very similar, but PSR 
J1227—4853 has a softer photon index and higher cutoff energy in this state than PSR J1023+0038. 

The gamma-ray light curve can be reproduced with outer-magnetospheric emission models for 
geometries that agree with the observation of radio eclipses and optical constraints on i, assuming 
that the orbital and spin axes have aligned through the recycling process. The near alignment of the 
gamma-ray peak with the main 1.4 GHz radio peak suggests that at least part of the radio emission 
may be generated in extended regions at high altitude in the magnetosphere, co-located with the 
gamma-ray emission site. With PSR J1824—2452A, this is the second gamma-ray MSP with high 
spin-down power (> 10 35 erg s _1 ) but only partially aligned light curve peaks, as opposed to PSRs 
J1823— 302 1A and B1937+21 for which the gamma-ray peaks are all aligned with radio features 
(e.g.. l,Tohnson et alJl2014l ). PSR B1957+20 represents a mix between the M SP subclass with aligned 
radio and gamma-ray peaks and the partially aligned MSPs. In particular, Guillemot et ah ( 20121 ) 
showed that at 0.35 GHz, both radio and gamma-ray peaks occur at the same rotational phase; 
however, an additional radio peak is present a t 1.4 GHz that is not matched by a feature in the 
gamma-ray light curve. Archibald et al. (120131 1 reported marginal evidence (3.7 a) for gamma-ray 
pulsations from PSR J1023+0038, with a suggested gamma-ray peak significantly separated in 
phase from the radio peaks (by 0.4 to 0.5 in phase). There is no evidence for modulation of the 
gamma-ray flux from PSR J1227—4853 at the orbital period during the pre- or post-transition 
periods. Continued study of the pulsed emission of PSR J1227—4853 such as increased gamma- 
ray statistics, radio polarization, and searches for pulsations at other wavelengths, should provide 
insights and help to further test and refine pulsar emission models. 

Coincident changes in the gamma-ray flux accompanied two of the three transitions in redback 
systems observed to date. In the case of PSR J1824— 24521, a n on-detection is not particularly 
troubling given the distance (5.1+0.5 kpc, Rees fc Cudworthlll991 .) and the added complication of 
a bright foreground signal in the LAT data from PSR J1824—2452A. As such, it is of interest to 
speculate on the rate at which redback systems undergo such transitions and whether we should 
expect to observe more transitions with the LAT. We know of nine redback systems in the Galactic 
field, which we can use to estimate a lower limit (acknowledging that there are unknown redbacks 
among the unassociated LAT sources) to the number of redback-years that LAT has observed to be 
54. With two events observed during the Fermi mission, we can estimate the rate at which a redback 
undergoes a transition to be once every > 27 years. We can check this estimate against the fact that 
PSR J1023+0038 showed evidence for an accretion disk in the 2000-2001 time fra me, was detected 


as a radio pulsar in 2009 , and transitioned back to an accreting state in 2013 ([Archibald et al 


2009 


Stappers et alJl2014J ). This suggests a timescale on the order of a decade between transitions, 
in rough agreement with our estimated lower limit. This lower limit on the rate may in fact 
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be overestimated if we consider that not all of the known redback systems in the Galactic held 
should be expected to transition at all. It is currently unclear what triggers a transition, in either 
direction. PSRs J1227—4853 and J1023+0038 both have P « 1.69 ms and orbital periods less than 
7 hours. If we only include redbacks with short spin (< 3 ms) and orbital periods (< 7 hours) 
we have four sources and our estimate becomes once every > 12 years, in very good agreement 
with the data from PSR J1023+0038. As the Fermi mission continues it is thus likely that other 
systems may transition, and monitoring the flux is important not only for known redback systems 
but also of unassociated sources with pulsar-like spectra from which pulsations have not yet been 
detected. Previous studies have defined pulsar-like sources to be non-variable with significantly 
curved spectra; however, given the measured spectral variability of PSR J1227—4853 during the 
LMXB phase it seems wise to relax the first criterion when searching for similar systems. 
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